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The role of shape resonances and many-body effects on universal quantum sticking of ultra cold 
atoms onto solid surfaces is examined analytically and computationally using an exactly solvable 
representation of the Dyson equation. We derive the self-energy renormalization of the the transition 
amplitude between an ultra cold scattering atom and the bound states on the surface in order to 
elucidate the role of virtual phonon exchanges in the limiting behavior of the sticking probability. We 

jr | ■ demonstrate that, to first order in the interactions for finite ranged atom-surface potentials, virtual 

phonons can only rescale the strength of the atom-surface coupling and do not rescale the range 
of the coupling. Thus, universal sticking behaviour at ultra-low energies is to be expected for all 
finite ranged potentials. We demonstrate that the onset of the universal sticking behavior depends 
greatly on the position of the shape resonance of the renormalized potential and for sufficiently low 
energy shape resonances, deviations from the universal s(E) oc y~E can occur near these energies. 
We believe that this accounts for many of the low energy sticking trends observed in the scattering 

^^ ' of sub-millikelvin H atoms from superfluid 4 He films. 



I. UNIVERSAL QUANTUM STICKING 



The exact limiting behavior of the quantum sticking probability has been repeatedly discussed over the past 25 
years, with various authors arriving at various conclusions |H,ppl|1.0,p|.Rj8l0] . Since the limiting behavior is expected 
only at extremely low scattering energies, only recently have direct experimental measurements of sticking at the 
required energies been obtained, most notably in the sticking of spin polarized hydrogen to liquid He films JlQ| , pT[ and 
in the desorption of cold positronium from metal surfaces [pL2| . At source of the theoretical controversy is the role 
O • that the many-body aspects of the scattering might play at extremely low energies and the effect of the long ranged 
^> behavior of the attractive part of the static potential. 

The first complete formal analysis of the problem was by Goodman in the early 1970's |jj. Using a WKB analysis 
of the low energy wavefunction, it was demonstrated that that lim^^o s(E) > for potentials with ranges longer 
than exponential and vanishes for those with ranges less than exponential even within the distorted wave Born 
approximation (DWBA) and that many-body effects are not needed in order to account for non- vanishing sticking 0] 
However, it must be pointed out that these results are semi-classical and that the WKB approximation used to derive 
their sticking results may not permit rigorous extrapolation of s(E) for energies below the WKB validity condition. 
Interestingly enough for attractive inverse power law potentials with a characteristic exponent r < 2, the WKB 
approximation is valid down to E = and classical and quantum predictions agree to within a constant factor. 0] 

In support of these results Knowles and Suhl p] have constructed a variational polaron model in which the scattering 
atom is dressed by its interaction with the surface phonons. Within such a basis, the adatom becomes a quasi-particle 
with an effective mass to* which is position dependent and much larger than the bare mass to over the potential region. 
This leads to an increase in the density of low energy resonances which increases penetration of the wavefunction 
into the potential region and leads to non-vanishing sticking even for square well potentials. Later, Brenig and co- 
workers y,|^,§| showed that within the DWBA: (1) sticking should always vanish for short ranged wells, (2) sticking 
can be enhanced by polarization induced resonances for longer ranged potentials, (3) for the Coulombic potential, 
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quantum and classical analysis yield identical results. Recently Carraro and Cole [J7[ applied Breni g's equations to 
study H sticking to liquid He film and obtained remarkable agreement with the experimental results fllf| . 

As far as we know, the most comprehensive examination of the quantum sticking problem is that of Cloughtery 
and Kohn who developed an analytic model for quantum sticking. [|S| Working up from a surface composed of N 
discrete lattice atoms to the continuum, N ~ * oo, they demonstrated that within their one dimensional model, while 
polarization effects due to virtual phonon exchanges may lead to an increase in the sticking coefficient, eventually 
these effects will be lost as pure reflection becomes more important at low energy. 

Experimental evidence for quantum sticking has come from a variety of different sources. Most notably in the 
desorption of positronium (Ps) from Al surfaces 12] and the sticking of spin polarized hydrogen atoms on liquid 4 He 
films |[(J,[L1|. Slow Ps atoms, composed of an electron and a positron, are almost ideal atoms for studying quantum 
sticking. At experimentally accessible scattering energies of 7.5 meV, the thermal de Broglie wavelength of Ps is on 
the order of 100 A and is much larger than the range of the potential well. Measurements of the thermal desorption 
rates of Ps from clean Al(lll) along with detailed balance arguments indicate that the Al surface is a "blackbody" 
for Ps emission and hence the system fails to exhibit perfect reflection of ultraslow Ps atoms. For the case of H on 
He films, the interactions are considerably weaker and there are a number of examples of experiments in which the 



sticking data extrapolates to s(0) = 1 1 3 14|jl5|,|l6|,|l7|,[l8 19 2flj21|. Surprisingly enough, recent data for H sticking 



onto thin He films, reported a gradual increase in the sticking probability as the temperature decreased between 10 
mK and .1 mK |10|Jll |, This prompted Hijmans, Walraven, and Shlyapnokov [B2j to propose that the atom-surface 
potential might be affected by very long-range van der Waals forces (with relativistic retardation) from the underlying 
substrate. Indeed, calculations by Hijmans, Walraven, and Shlyapnokov p2| and by Carraro and Cole H indicate 
that substrate effects could account for this trend and predict that the y/T behavior will be recovered at slightly lower 
temperatures. More recent work which measured the sticking coefficient as a function of the He film thickness gave 
reasonable agreement with the theoretical predictions in temperature regimes above .1 mK [[nj That is, as the film 
thickness was increased, the sticking probability was reduced and eventually the universal s(T) oc VT behavior was 
observed. 

In this paper we wish to examine two issues which we believe have not been adequately addressed, and yet are 
necessary in order to properly interpret the experimental results mentioned above. The first issue is the role of 
many-body effects on the limiting behavior, and the second is the role of low energy shape resonances. In the 
next section, we examine the many-body contributions by solving formally the Lippmann-Schwinger equation for 
the scattering wavefunction in which we include the self-energy due to the virtual phonon transitions between the 
scattering wave and the inelastic channels. Since our theory incorporates the dynamical evolution of the surface 
directly into the final equations of motion, it is essentially exact (for single phonon exchanges) and non-perturbative. 
Using approximate forms of the low energy scattering wavefunction, which we show to become exact as k — ► 0, we 
derive the renormalization of the sticking probability due to the the many-body interactions for realistic atom-surface 
interactions. We then demonstrate that polarization due to many-body effects serves only to rescale the strength of 
the interactions and does not rescale the effective range. We next examine the effect of low lying shape resonances and 
demonstrate the effect of increased penetration of the scattering wave into the potential region near these resonances. 
Finally, in the last section we apply the methods which we developed to examine the low energy limiting behavior of 
H sticking to thin and thick 4 He films. 

II. THEORY 

The starting point of our theory is the Lippmann-Schwinger (or Dyson) equation for the scattering wavefunction in 
which we include explicitly the self-energy or polarization due to virtual phonon exchanges. The self-energy reflects 
the many-body nature of the problem and we shall assume that only single exchange processes are important. Close 
coupled forms of this equation have been used extensively by Stiles, et. al p3|,|24| and Jackson [ p5p^j27| , |28| , p9| , p0] , pT| to 
study inelastic molecule-surface scattering and Whaley and Bennett have used the formalism to study atom scattering 
from disordered surfaces |^1 ■ We write the elastically scattered component of the wavefunction as 

ME) = 4>bare{E) + G bare (E)V sel f(E,T)4> {E), (1) 

where <fibare(E) and Gbare are the bare wavefunction and Green's functions in the absence of surface vibrations and 
Vseif (E, T) is the self energy operator 

V self (E,T) = ^Y, V t ( G+ q (E,T) +G-(E,T))V? h . (2) 



Our notation is such that superscripted + and — refer to whether or not a virtual phonon was initially created or 
annihilated. Also, throughout this paper we shall define H as the static (zero-phonon) Hamiltonian and V^ h (z) 
as the force between the scattering atom and a phonon with wavevector q. We have also assumed that the bath is 
harmonic. 

Closure is accomplished in the usual way by iteratively substituting <f> back into the Lippmann-Schwinger equation 
and analytically performing the summation over single phonon exchange diagrams. 

4>o = <Pbare + GbareV S elf4>bare H 

= (1 — GbareV se lf) §bare (3) 

Before applying our theory to a physically realistic system, it is important to explore the effects induced by the 
additional polarization due to virtual phonon exchanges. Let us consider the "golden rule" transition rate between 
a very low energy scattering state and a single bound state, ip B , and attempt to relate the transition probabilities 
predicted using a dressed wave to those predicted from using the bare wavefunction. In short, we need to consider 
the matrix elements 

(MV h \4>o) = {M {1 - V ph G {E)V^ h {G + {E) + G-iE))}- 1 ^|</w), (4) 

where G (E) is a solution of 

(H -E)G {z,z')=6{z-z') (5) 

subject to the boundary condition G (— oo, z) — and G (z, z') — > as z — > oo. 

To evaluate this matrix element we need to make a series of approximations regarding the functional form of 
4>bare- The bare scattering wave, 4>bare, is an solution of (H — E)<fibare — subject to the boundary conditions: 
4>bare ( — oo) = and <fibare (z) — + sin(kz) + tan(<5fe) cos(fc;z) as z — > +00. We can also define the real regular and 
irregular functions, ip r ,i(z), as solutions of 



(H o -E)il> r 4z;E) = (6) 



subject to the boundary conditions 



ip r (z,k) — ► sin(fcz + S) 
ipi(z,k) — ► cos(kz + 5). 

(7) 

Using these functions and the Wronskien relation, we can write G (z, z'\ k) as 

r (■? yi-h\- _i / ^r&fyi/jiiz'ik) +iipr(z;k)ipi(z';k) z < z' , ■> 



where k = y2mE is the scattering momentum. 

For very low energies and temperatures, virtual transitions are predominantly to the bound state and we can 
obtain the single phonon exchange propagator from the bound (asymptotically closed) wavefunctions, weighted by 
the appropriate phonon densities of states, which in the continuum limit can be taken as 

G+(z, z>) = M z) M z') I p. ^gCgg/q (Wg (D + 1) 

J 2tt uj{q){E B -u{q) + k 2 ) 

= M'WbW r dupiu) C f {q ^ )a W ( W ( W , T) + 1). (9) 

J \k z + E B - U) + iri)LU 

where lud is the Debye frequency of the lattice, n q (T) and n(u,T) are the Bose-Einstein occupation numbers, and 
p(u>) is the density of states. The integral can be evaluated by choosing the appropriate contour which avoids the 
pole at uj = k 2 + E B ra and the result is (using a Debye model for the phonons) 

G+(z,z')=^ B (z)^ B (z')J'(k), (10) 

where 
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(11) 



The notation -f, dx means that we take the Cauchy principal value when evaluating the integral. 
Putting all this together yields an expression for the matrix element which is convenient to evaluate 



dz / dz'^ B {z){z\l-V ph G V ph Gi\z')- l V ph {z')<j> bare {z'). 



(Mv ph \^o) ■- 

Here, self energy matrix element is given by 

(z\l - V ph G V ph G+\z') = S(z - z') - V ph (z)i> B (z')T(k)I{z- fc), 
where 



I(z;k) 



dz"G {z lZ "-k)V ph {z")^B{z") 



(12) 



(13) 



(14) 



In order to take the k — > limit of these expressions, we first need to make an assumption regarding the extreme 
low energy behavior of the ipo t i(z; fc) functions used to construct G . (Justification for this assumption for the 1/ z 3 
potential is included in the Appendix.) 



lim ipr(z~, fc) ~ ktjj r (z, k = 0) 
lim ipi(z; fc) w ipi(z, fc = 0) 

fe— »-0 



(15) 
(16) 



where ip rt i(z, k = 0) are solutions of H ipQ^i(z) = which are regular and irregular at the surface and are independent 
of fc. Using these, we can construct the E — > Green's function 



lim G (z. z'\ fc) 

fc^O 



ij} r {z; k = Q)tpi(z'; k = 0) + ikip r (z; k = 0)ip r (z'; fc = 0) z < z' 
i> r (z'; k = 0)t/ji(z; k = 0) + ikip r (z; k = 0)ip r (z'; k = 0) z > z' 



(17) 



Thus the integral, I(z; fc), which occurs in the self-energy expression can be broken into real and imaginary components, 
both independent of fc. That is 



(z|l - V ph G V ph G+\z') h =° 5{z - z') - V ph {z)^ B {z'){h{z) + ikl 2 {z))f(k), 



where 



(18) 



(19) 

(20) 



(21) 



(22) 

Since the integral in £(fc) is at most linearly dependent on fc, any deviations from the DWBA results are likely to 
come from T{k). We note that the denominator, 1 — S(fc), should approach a non-zero limit as fc —> 0. The limiting 
low energy behaviour will thus be determined by the perturbation integral, (ifjB\V ph \(j)bare), which approaches zero as 
\/E in the limit of fc — > for finite ranged atom-phonon interactions. 

What we have neglected in this evaluation is the possible contribution from shape and threshold resonances. In 
particular, a threshold resonance or a shape resonance nearly at the threshold may produce an additional pole close 



h(z) = ftel(z:k^0) 
h{z) = %mI(z- 1 k->0). 

Finally, substituting this back into Eq. H yields for small fc 

(tPbIV^I^o) = (yj B \V ph \fo are )(l + E(fc) + S 2 (fc) + • • •) 

= (lpB\V ph \4> ba re)/(l-^(k)), 

where £(fc) is the polarization renormalization due to virtual phonon exchanges ( in the small fc limit). 

E(*)=^(Jb) f dzV ph (z)(I 1 (z) + ikI 2 {z)). 



enough to the real energy axis to change the limiting behavior of T(k). For the moment, let us consider the case that 
there are no threshold resonances and that we are well below the last shape resonance of the static potential surface. 
In this case, E(0) is a constant between 1 and -1 and Eq. |2l] produces only a rescaling correction to the DWBA result, 
i.e. there is no change in the effective range of the potential and the s — ► yE limiting behaviour should hold. 

Interestingly enough, the calculations which we presented in Ref. H , did not show the limiting behavior predicted 
above and we have identified the reasons. In that paper, we computed the sticking probability as a function of surface 
temperature at a scattering energy E — > for a model He atom scattering from a featureless Cu surface. We first 
tested the model by comparing finite E results to results published previously for the same model and obtained good 
agreement. We then varied both the attractive range of the potential, from long ranged V a tt oc — l/z 3 to the short 
ranged V a tt °c —e~ az and the atom-phonon interaction by varying the phonon force constants and by changing the 
effective mass of the surface atoms. We then reported that the model produces sticking as E — ► which decreases as 
the range of the potential is decreased and decreases as the phonon force constants are increased. Since this is precisely 
the behavior one would expect for potentials with longer ranged interactions or if the effective scattering mass of the 
atom is increased, we concluded that the inclusion of the virtual phonon processes causes a renormalization of the 
range of the static potential and/or an enhancement of the effective mass of the scattering atom and, hence, leads to 
non- vanishing sticking. 

Later analysis indicates that these calculations may be in error due to our particular choice of a linear trial function 
used in our numerical calculations to generate the E — ► solutions. In our original calculations, we tested the method 
against results published by Jackson for He scattering from Cu(100) at 5.1 meV and normal incidence and found 
satisfactory agreement and convergence. The problem with the low E results occurs because the finite ranged basis 
and the linear auxiliary trial function used imposes an artificial energy minimum in the convergence properties of the 
equations and the method can not converge wavefunctions with energies below that critical value. Usually, convergence 
errors occur at the high end of the energy spectrum where one does not have enough high energy components in the 
basis used to represent the Green's function and thus to adequately represent the wavefunction. In the present 
calculations, however, we use sin(fcz) as the auxiliary trial function and do obtain the correct limiting behaviors 
predicted by the above analytical model. 

In the next section, we apply the equations derived above to study H reflection from He films. Since hydrogen is 
an extremely light atom and the He-H potential is well characterized this is an ideal physical system for this type of 
study. Also, it has been demonstrated recently, that the sticking behavior is strongly dependent upon the long ranged 
interactions with the underlying substrate p"l|M|22|]. We shall demonstrate that the increased polarization from the 
underlying substrate causes an increase in the density of low lying shape resonances which enhance penetration of the 
wavefunction into the potential region and enhance sticking. To elucidate the effect of resonances in the sticking, we 
examine sticking on an even more idealized situation in which the H-He interaction is reduced to that of a square well 
and H-ripplon interactions decay exponentially in z. By adjusting the "range" of the well (at constant well depth) 
we can very easily see the enhancement of sticking due to low lying resonances. In contrast to our earlier results, 
once the scattering energy is sufficiently below the resonances, penetration into the potential region becomes greatly 
diminished as reflection from the long ranged part of the potential begins to dominate and the universal s(E) oc y/E 
behavior of the sticking coefficient is recovered. 

III. STICKING OF ATOMIC HYDROGEN ON LIQUID HELIUM FILMS 

The interactions between H and liquid He are known rather well. In fact, variational estimates of the binding 
energy of H on liquid 3 He and 4 He date back to the late 70's |I^,|l3[. At long range, the potential can be obtained 
from summing over all pairwise van der Waals interactions between the H and the He atoms. At large distances 
(over 200A), retardation effects need to be included in order to provide the proper cut off in the attractive forces. At 
shorter ranges, the potential can obtained from empirical fits to atomic beam scattering data. |lSfl . To lowest order, 
the effective H-liquid He potential can be constructed by integrating over H-He pairwise potentials and assuming that 
the liquid He is a semi-infinite and incompressible fluid extending from z = — oo to z = and with a sharp density 
profile, p(R, z) — p Q{z + u(R)). 

r r + co 

V(z) = Po / d 2 R' / dz'v pair (r)Q(z' + R'), (23) 



where r = \J R 2 + (z — z') 2 , v pa i r is the H-He pair potential, and p is the number density of bulk liquid 4 He. 
( / 9 = 0.0218A- 1 .) 

There have been numerous theoretical treatments of the H-He pair potential, most recently by Bhattacharya and 
Anderson using a quantum Monte Carlo method and the reader is referred to their paper for a comprehensive listing 



of previous treatments [[33[ . A number of forms for the effective potential have been proposed ranging from Morse 
potentials used by Zimmerman and Berlinski |13] to much more realistic forms which include the 1/z 3 van der Waals 
contributions both from the He film and the underlying substrate. The effective potential used in our study is similar 
to that used by Carraro and Cole H and by others |p4| . We included both He film and substrate van der Waals 
polarization contributions as well as a relativistic cut off factor which limits the long range behavior of the potential 
at distances greater than A = 200A. 

v ° = A(z) - Cs ( (*-*!)«+*! - jdw) 7(z) " Cs idw l{z + d) ' (24) 

where A(z) = Dexp(— (3{z — z )) and arises from atomic core repulsions. As discussed in the introduction, the final 
results should not be very sensitive to the exact form of this part of the potential. Retardation effects in the attractive 
terms are included in 7(2) which we assume to take the form 7(2) = (1 + z/A) -1 where A = 200A is the characteristic 
range of the effect. The remaining parameters, 0, z and D are tabulated in Table | and were constrained to produce 
a single bound state with a binding energy within the range of the latest experimental estimates of 1.0 ± .1 K for an 
infinitely thick He film. The polarization terms are C3 = 219.7 K A 3 and C s = 5000 K A 3 . Finally, we note that any 
effective potential, V, can be generated by an effective two-body potential, v by directly inverting Eq. G3. pM 

We treated the He film as a smooth structureless surface. In this system, the inelastic effects are mediated through 
the capillary motion of the film. For the temperature ranges in which we are interested, these capillary waves, called 
ripplons, provide the only remaining inelastic channels, the other inelastic channels, such as those mediated by the 
phonons and the rotons are effectively frozen out. Mathematically, the ripplons and phonons are essentially equivalent 
quantities, differing only in dispersion and coupling. The displacement of the He surface can be written in terms of 
the ripplon operators as 

MK> = ^£(^) 1,2 <4 + «-,>^. (26) 

The ripplon dispersion for thick He films is given by the well know form for surface capillary waves 

^ = ^ + w P q3 ' (27) 

where a is the He surface tension (a = 0.27 KA~ 2 ) pa], and the "gravitational" acceleration is given by, g = C s /d A . 
As with He scattering from metal surfaces, we make the ansatz that the complicated dependence upon q in the 
atom-ripplon interaction can be factored into two terms, one which includes all information regarding momentum 
transfer and one which depends upon only the core interactions at the surface. The form which we use is 



The fact that we are using an approximate form for the coupling may introduce an overall error in our final results; 
however, the overall limiting behavior (i.e. whether or not the sticking vanishes) is expected to be universal and 
should not depend strongly upon the exact form of this coupling. 

Finally, the sticking probability is obtained from the transition probability between the dressed scattering state and 
the bound states. 

S (T,E) = ^ J2 l<W l^>l VW& - w«)/(l " e"^) (29) 

n 

where if)" (uj) is the imaginary part of the ripplon susceptibility, which can be approximated using a Debye model, 

1>'>) = 2j£jSve(vD-\w\), (30) 

where u>d is the surface Debye frequency, which we take to be on the order of 1 meV, or about 12 K. This produces 
a cut off in the ripplon momentum dispersion of qo = 1 A -1 . 



Figure |l| shows the variation of s(T) for a 30 A He film and bulk He as the H scattering energy is reduced below 10 
mK. Superimposed are the experimental data for the same temperature regime. Given the uncertainties in both the 
experimental data and in the sensitivity of our calculation to variations in the potential parameters, the agreement is 
somewhat satisfactory. For scattering from bulk He films, the substrate polarization is attenuated by the thickness of 
the film. Although the binding energy of H on bulk He is somewhat less than that for thin films, the shape resonances 
are far enough into the continuum that penetration into the potential region is not enhanced below 1 mK and quantum 
reflection begins to dominate. For thin films, the situation is just the opposite. The additional polarization due to 
the underlying substrate dominates the long ranged behavior of the potential and the density of the shape resonances 
near threshold is increased. This gives rise to the observed increase in the sticking as the film thickness is decreased. 
Below about 0.1 mK, which is at the end of the most recent available sticking data, barring the presence of a lower 
lying shape resonance, the resonance enhancement effect is diminished and quantum reflection should be observed. 

In the next figure, (Figure ||), we show the changes in s(T) as the He film thickness is increased from thin film to 
bulk at a constant scattering energy of .35 mK. Again, although our results differ from the experimental values by 
perhaps a factor of two or three on the average, the agreement between our calculations and the experimental data 
is on the whole quite good, especially given the apparent sensitivity of s(T) to the long ranged part of the potential 
and the lack of an accurate measurement of the substrate interactions. What is striking about our calculations is the 
bump in our computed s(T) which appears when the film thickness is just below 50 A. This is due to the appearance 
of a threshold resonance, which eventually changes to a second bound state when the film thickness is around 15- 
17 A. In Table 1 we show the variation in the binding energy as the film thickness is decreased as computed by 
diagonalizing the static surface Hamiltonian. Quantitative comparison to the the experimental data much below this 
region is difficult due to other substrate effects which might come into play, such as roughening, corrugation effects, 
and inelastic coupling to the substrate phonon modes. In fact, the experiments were done using either sintered silver 
or epoxy as a substrate. For the sintered silver case, for coverages above 30 A, the pores in the sinter should be 
filled jy] ; however, for coverages below that the pores may be only partially filled. In any case, the experimental data 
also indicates an abrupt change in the sticking behavior which occurs when the He film thickness is below 20 to 30A. 

As a final demonstration of the role of shape resonances in quantum sticking, we consider the case of H sticking 
to "pseudo" liquid He where we have replaced the realistic effective potential used above with a square well. By 
adjusting the width of this well, A, we can very easily examine (analytically, if we want), the enhancements to the 
sticking due to low lying resonances. In Fig. @ we plot s(T) vs A at constant scattering energy well into the quantum 
reflection regime (E — 0.01 meV). For the narrowest well (A = 10 A), the potential barely supports a single bound 
state with a binding energy of 0.031558 K. As A is increased to 50A, more bound states are added to the well. This 
affects the sticking in a profound way. At each peak in this figure, an additional bound state is added to the well. As 
A increases further, the binding energy of the new state increases and sticking decreases. 

Lastly, in Fig. [l| we consider a rather exceptional case in which the effective range of the atom-phonon interaction 
was chosen to be longer than the effective range of the attractive part of the potential. Here, we set 1/ '(3 = 20A 
and vary both the range of the square well, A, and the scattering energy. Two features are immediatly appearent. 
First is that the universal sticking behaviour is observed in all cases indicating that for finite ranged atom-phonon 
interactions, universal sticking is to be expected. Secondly, is the dramatic enhancement of the sticking due to the 
low lying resonances. 

IV. DISCUSSION 

We have presented a theory for the sticking of ultra-cold atoms from surfaces, and we have presented evidence 
for what we believe to be the proper limiting behavior of s(T) for non-Coulombic potentials. Although our model 
calculations are limited to low scattering energy where only transitions to bound states are allowed, we believe that the 
limiting behavior shown here is in fact correct, agreeing with the predicted behaviors. Our calculations indicate that 
s(T) can have non- vanishing limiting behavior only in the presence of extremely low lying (threshold) resonances. 
We demonstrate that many body effects do persist into the quantum sticking regime; however, their net effect at 
very low energy is to rescale the strength of the potential but not to rescale the effective range. Since sticking is 
most sensitive to the range of the potential, the limiting low E behavior becomes independent of the strength of the 
interaction and the coupling to the inelastic channels. At very low scattering energy, there is a competition between 
polarization effects due to the virtual phonon exchanges and quantum mechanical reflection by the long ranged part 
of the potential. In the absence of shape and threshold resonances, quantum reflection is expected to win out and the 
universal s(E) ex \/E limit is obtained. 

There are a number of fundamental questions which remain unresolved both experimentally as well as theoretically. 
When the inelastic interactions are very weak, as in the case of H on He films, subtle details, such as relativistic 



retardation effects and substrate contributions to the potential become very important and can dramatically change 
the sticking behavior. These subtle contributions are very difficult to quantify experimentally due in part to the 
technical difficulties of working in the submillikelvin regime [Q. However, given the sensitivity of the sticking to very 
subtle details of the long range part of the potential, these low energy experiments offer unique insight to substrate 
effects and relativistic retardation effects that would otherwise be unmeasurable. 
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APPENDIX A: ZERO ENERGY WAVEFUNCTIONS FOR 1/Z 3 POTENTIALS 

The essence of the theory of quantum sticking at very low scattering energies boils down to the fact that near the 
classical turning point, the scattering wavefunction vanishes linearly with momentum, k. Thus, any matrix element 
coupling this low energy scattering state to a bound state will also vanish. A demonstration of this can be seen by 
solving the following Schrodinger equation, written in dimensionalless units. 



1% + (g z /z A + k z )^ k = 



(Al) 



on z > 1 and subject to the the boundary condition, ipk(X) = 0. Exact analytic solutions for this exist only for k = 
and the result is [5[ 



M*) = Ay/z[Ni(x) ~ N 1 (2g)/J 1 (2g)J 1 (x)} , 



(A2) 



where Ni(x) and J\(x) are Neumann and Bessel functions with arguments, x = 2g/w(z). The asymptotic expansion 
of rpo is obtained by taking x — > and using the known asymptotic expansions for the Bessel and Neumann functions 



J x ( x ) = -- — 

V ' 2 16 



Ni(x) 



xln( , yx/2) 



(A3) 

(A4) 



where 7 is the Euler constant (In 7 = .5772 • ■ ■). Also, we can use the x — > 00 asymptotic forms of J\ and N\ to obtain 



N l (2g)/J 1 (2g) = tan(2g- -tt). 



(A5) 



Putting things together, one obtains 

i>o(z) = ~Ag 



tan< 



In 



■ ! \1 2 9 2 



,9 Z 



(A6) 



For k y^ 0, approximate solutions can be obtained from the k — solution by partitioning the potential into interior 
and exterior portions. The interior part, 



Vo(z) 



g 2 1 ' z 3, — k 2 forz < z c 
otherwise, 



where the critical distance, z c is determined by k 2 = g / z 3 . From this particular choice, it is obvious that for z < z c , 
the above expression for tpo is the solution, and for z > z c , the asymptotic form, 



ipk = sin(kz c + 5 k ), 



(A7) 



is the solution. The phase shift, Sk, and coefficient, A, are determined by requiring the the logderivative, 7fc(z) = V'fc/V'fc 
be continuous at z = z c . From the matching condition, as k —> 0, 



A = -7T gk + 0[k 5/2 Ink}. (A8) 

Similarly the phase shift can be estimated. The important observation is that the amplitude of the wavefunction inside 
the interaction region, as defined by z < z c , vanishes linearly with k. This is to say, that there is no transmission of 
probability flux past z c , and thus the wavefunction is totally reflected by the attractive part of the potential before it 
even reaches the spatial region in which it can interact strongly enough with the phonons to cause sticking to occur. 
This is the origin of universal quantum reflection. 
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TABLE I. Effective Potential Parameters and Computed Binding Energies for H on Bulk 4 He and on thin 4 He films (d = 
30 A). The reported value of E„ for bulk 4 He is 1.0K (See Ref. jnj.) 

D[K] (5 [A' 1 ] Zo [A] E [K](thin) E [K] (bulk) 

3.0 0.20 4.8 -0.981589 -0.855165 

2.9 0.20 4.8 -1.059139 -0.931325 

2.8 0.19 4.8 -1.155515 -1.025170 

2.8 0.20 4.8 -1.139079 -1.009945 

2.8 0.21 4.8 -1.122955 -0.994995 

2.7 0.20 4.8 -1.221373 -1.090980 

2.6 0.20 4.8 -1.305987 -1.174390 



FIG. 1. Trapping probability versus atom energy for H onto a 3 nm thick and bulk He films (1 mm). The curves labeled A 
and B are the results obtained from our calculations for thin and thick He films. Superimposed are the experimental results 
for each of these limits. 

FIG. 2. Sticking probability versus He film thickness. In each case, the scattering energy and surface temperature was .35 
mK. The diamonds are the experimental data points from Ref. [ |Q0j] The bump in our computed data is due to a threshold 
resonance which appears when the film is less than about 30 A thick. 

FIG. 3. Sticking vs. well width for a square well potential at constant scattering energy, (well depth = 2.8 K) The increases 
in the sticking are due to the the appearence of additional weakly bound states as the well width is increased. For the narrowest 
case considered, the well supports a single bound state almost at threshold (binding energy = 0.03 K). Threshold well widths 
for each additional bound state are indicated by arrows. 

FIG. 4. Sticking vs. Scattering energy for a square well potential. Here we have chosen the effective range of the atom-phonon 
to be 20A and vary the width parameter, A from 9.5A to 21. 5A in which we pass from 1 bound state to 2 bound states in the 
square well potential. See text for details. 
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